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Abstract 

We consider a scalar dark matter model, the SM4+D, consisting of the standard model with four 
generations (SM4) and a real gauge-singlet scalar called darkon, D, as the weakly interacting massive 
particle (WIMP) dark-matter (DM) candidate. We explore constraints on the darkon sector of the 
SM4+D from WIMP DM direct-search experiments, including CDMS II and CoGeNT, and from the 
decay of a B meson into a kaon plus missing energy. We find that a sizable portion of the darkon 
parameter space is still compatible with the experimental data. Since the darkon-Higgs interaction 
may give rise to considerable enhancement of the Higgs invisible decay mode, the existence of the 
darkon could lead to the weakening or evasion of some of the restrictions on the Higgs mass in the 
presence of fourth-generation quarks. In addition, it can affect the flavor-changing decays of these new 
heavy quarks into a lighter quark and the Higgs boson, as the Higgs may subsequently decay invisibly. 
Therefore we also study these flavor-changing neutral transitions involving the darkon, as well as the 
corresponding top-quark decay t — > cDD, some of which may be observable at the Tevatron or LHC 
and thus provide additional tests for the SM4+D. 
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I. INTRODUCTION 



The existence of dark matter (DM) in the Universe is now widely accepted. Various observa- 
tions have established that DM makes up 23% of the total cosmic energy density [ij. Despite this 
evidence, however, the identity of the basic constituents of DM has so far remained a mystery. 
It is therefore important to explore different possible DM scenarios. 

One of the popular candidates for DM is the weakly interacting massive particle (WIMP). 
To account for WIMP DM, the standard model (SM) of particle physics needs to be enlarged. 
The simplest extension of the SM possessing a WIMP candidate is the SM+D, which combines 
the SM with a real SM-singlet scalar field D, dubbed darkon, to play the role of the DM. This 
darkon model and some variations of it have been much studied in the literature (2-12 . 



In this paper we explore a somewhat enlarged darkon model we call SM4+D, which con- 
sists of the darkon and the SM extended by the inclusion of a fourth sequential generation of 
quarks and leptons. This SM with four generations (SM4) has received much attention in recent 



years [13H23J. Among the reasons [13| that have been put forward for all this interest in the 
SM4 are that it is not ruled out by electroweak precision tests 14j-|l6l|. offers possible resolutions 
for certain anomalies in flavor-changing processes 17H20||. and might solve baryogenesis-related 



problems [21]. In view of the desirable features of the model, some of which nevertheless remain 
open questions, it is of interest also to consider integrating the darkon field into it, assuming 
that the new fermions are all unstable, in which case the SM4+D is the simplest WIMP DM 
model in the presence of the fourth generation. As we will elaborate later, the DM sector of 
the SM4+D can have important implications which are absent or suppressed in the SM+D with 
three generations (hereafter referred to as SM3+D). In particular, now that the LHC is opera- 
tional, the extra fermions could give rise to processes involving the darkon which are potentially 
observable after the LHC reaches full capacity in the near future. 

In the next section we describe the main features of the SM4+D relevant to our study. 
Subsequently, after specifying the masses of the fourth-generation fermions, we extract the values 
of the darkon-Higgs coupling, to be used in later sections. In Sec. lIIII we explore constraints on 
this darkon model from DM direct searches at underground facilities. Recently there have been 
a number of such searches which can provide limits on some of the parameter space of the 
darkon model. We proceed in Sec. lIVI to discuss the complementarity of DM direct-detection 
experiments and Higgs studies at colliders in probing the darkon properties. The simultaneous 
existence of the darkon and 4th-generation fermions in the SM4+D can have substantial impact 
on Higgs collider searches. Since ongoing and near-future DM direct-search experiments are not 
likely to be sensitive to darkon masses of a few GeV or less, other processes are needed to probe 
the model in this low-mass region. In Sec. [V] we consider such processes, focusing on the 5-meson 
decay into a kaon and a pair of darkons, B — > KDD, which contributes to the B decay into K 
plus missing energy, B — > K$. There is currently experimental information on the latter decay 
which can be used to place restrictions on part of the darkon low-mass region. In Sec. lVIl we 
explore some implications of the new fermions for the darkon sector that are lacking or missing 
in the SM3+D. Specifically, we look at the Higgs-mediated flavor-changing top-quark decay 
t — > cDD, which is very suppressed in the SM3+D and can be greatly enhanced by the new- 
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quark contribution, and also deal with the corresponding decays of the 4th-generation quarks. 
These processes may be detectable at currently running or future colliders and, if observed, could 
offer additional means to probe darkon masses from zero up to hundreds of GeV. We give our 
conclusions in Sec. lVIII 

Before proceeding to the next section, we would like to summarize the relic-density require- 
ments that any WIMP candidate has to meet. For a given interaction of the WIMP with SM4 
particles, its annihilation rate into the latter and its relic density Q D can be calculated and are 
related to each other by the thermal dynamics of the Universe within the standard big-bang 
cosmology jiij . To a good approximation, 



„ 1.07xl0 9 x f 0.038 m P1 m n (a v ,) 



where h is the Hubble constant in units of 100 km/(s-Mpc), m P1 = 1.22 x 10 19 GeV is the Planck 
mass, Xj = m D /Tj with Tj- being the freezing temperature, is the number of relativistic 
degrees of freedom with masses less than Tj, and (o" ann f rcl ) is the thermally averaged product 
of the annihilation cross-section of a pair of WIMPs into SM4 particles and the relative speed of 
the WIMP pair in their center-of-mass (cm) frame. Since Qd is known from observations, using 
the above relations one can extract the allowed range of <r ann for each value of m D . 



II. BRIEF DESCRIPTION OF SM4+D 

Being a WIMP DM candidate, the darkon D has to be stable against decaying into SM4 
particles. This can be realized by assuming D to be a singlet under the SM4 gauge groups and 
introducing a discrete Z 2 symmetry into the model. Under the Z 2 transformation, D — > —D, 
while all SM4 fields are unchanged. Requiring, in addition, that the darkon interactions be 
renormalizable implies that D can interact with the SM4 fields only through its coupling to the 
Higgs-doublet field H . It follows that the general form of the darkon Lagrangian, besides the 
kinetic part \d^D d„D and the SM4 terms, can be expressed as [3-0] 

C D = -^D*-^D 2 -\D 2 H^H , (2) 
4 2 K J 



where X D , m , and A are free parameters, and we have followed the notation of Ref. 10j. The 
parameters in the potential should be chosen such that D does not develop a vacuum expectation 
value and the Z 2 symmetry is not broken, which will ensure that the darkon does not mix with 
the Higgs field, avoiding possible fast decays into other SM4 particles. 

The Lagrangian in Eq. (J2]) can be rewritten to describe the interaction of the physical Higgs 
boson h with the darkon as 1 

C D = - ?f D* - (m ° + XV ^ D 2 -±D 2 h 2 -\vD 2 h, (3) 



1 Obviously, h here is not to be confused with the Hubble constant, also denoted by h, in the combination iljyh 2 . 
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where v = 246 GeV is the vacuum expectation value of H. The second term in C D contains 
the darkon mass m D = fm| + Aw 2 ) 1 / 2 , and the last term, —XvD 2 h, has a major role in the 
determination of relic density of the darkon. Clearly this model has a small number of free 
parameters in its DM sector: the darkon mass m D , the Higgs-darkon coupling A, and the darkon 
self-interaction coupling A D , besides the Higgs mass m h . Our analysis will not involve X D . 

For m D < m h the relic density results, at leading order, from the annihilation of a darkon 
pair into SM4 particles via Higgs exchange j^-ji]], namely DD — > h* — > X, where X indicates 
SM4 particles. Since the darkon is cold DM, its speed is nonrelativistic, and so a darkon pair 
has an invariant mass \/s ~ 2m D . With the SM4+D Lagrangian determined^the h- mediated 
annihilation cross-section of a darkon pair into SM4 particles is then given by [4| 

8AV E,r(ft->*) 

(iml-mlf + rlml 2m D ■ U 

where v rel = 2\pJ^ l \/m D is the relative speed of the DD pair in their cm frame, h is a virtual 
Higgs boson having the same couplings to other states as the physical h of mass m h , but with 
an invariant mass y/s = 2m D , and h — > Xj is any kinematically possible decay mode of h. 
To determine T^T^h — > XA, one computes the h width and then sets m h equal to 2m D . For 
m D > m h , darkon annihilation into a pair of Higgs bosons, DD — > hh, also contributes to o" ann , 
through s-, t-, u-channel, and contact diagrams with vertices arising from the last two terms of 
Cd in Eq. ([3]) and the Higgs self- interact ion (see, e.g., Ref. (js)]). This becomes one of the leading 
contributions to a ann , along with DD -)■ h* — > WW, ZZ, if m D > m WZ h [!, 0]. 

Compared to the SM3+D case, one effect of the fourth generation of quarks and leptons in 
the SM4+D is to enlarge the Higgs total width, T h , and also the total width HiT(h — > XA of 
the virtual Higgs, h. These new heavy fermions contribute to the total widths mainly via the 
decay modes into fermion-antifermion pairs if kinematically possible and, exclusively for the new 
quarks, the decay mode into a gluon pair induced by a quark loop. Needless to say, the changes 
caused by the presence of these fermions depend on their masses. 

There are constraints on the masses of the 4th-generation fermions from currently available 
experimental data. The masses of the heavy charged-lepton £' and heavy neutrino z/, both 
assumed to be unstable, have the PDG lower bounds m i , > 100.8 GeV and m v , > 90.3 GeV 
For the masses of the new up- and down-type quarks, if and b', respectively, the strongest limits 



are m t , > 311 GeV and m b , > 338 GeV, from searches at the Tevatron [251 ] . The mass differences 
between the new quarks and between the new leptons turn out to be subject to empirical 
constraints as well. Electroweak precision data prefer m t ,—m b , ~ [5+ln(m fo /115 GeV)] x 10 GeV 
and 30 GeV < m f , — m v , < 60 GeV [3] ■ Accordingly, for numerical work in this paper we take 
for defmiteness m e , = 200 GeV, m v , = 150 GeV, m t , = m b , + 55 GeV, and m t , = 500 GeV, but 
we also sometimes make comparisons with the m t , = 400 and 600 GeV cases. We remark that 
these m t , values fall within the ranges allowed by recent global fits for the SM4 15Ml9||. although 



m t i = 600 GeV is slightly above the unitarity upper-bound of ~550 GeV [13[ . 

With these mass choices, we can find the Higgs total widths, which we subsequently apply in 
Eq. (j3J), combined with the DD — > hh contribution if m D > m h , in order to extract the darkon- 
Higgs coupling A for given values of m D , m h , and (cr^v^). The allowed range of (cr ann t> rol ) as 
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a function of m D can be inferred, with the aid of Eq. (JTJ), from the data on the relic density. Its 
most recent value is Q D h 2 = 0.1123 ± 0.0035, determined by an analysis of the seven-year data 



from WMAP combined with other data [26J. From this number, one can derive the 90%-C.L. 
range 0.1065 < Q D h 2 < 0.1181, which we adopt for our numerical study. We show in Fig.[TJ 
the resulting ranges of A, taken to be positive, corresponding to 3 GeV < m D < 1 TeV for 
some specific values of the Higgs mass, which we choose to be m h = 115, 200, and 300 GeV 
for illustration. We present plots both in the SM3+D and in the SM4+D with m t , = 500 GeV 
for comparison purposes. The SM4+D plots for m t , = 400 and 600 GeV turn out to be very 
similar to the one displayed. 

There are several points worth pointing out in relation to what can be seen in Fig.[TJ First, 
although only a relatively narrow range of the DM relic density is allowed, evidently it can be 
fairly easily reproduced in both the SM3+D and SM4+D. Second, A is not small for the lower 
values of m D , and this will result in a considerable branching ratio of the Higgs invisible decay 
mode in the two models, as we will discuss further later. Third, for 3 GeV < m D < 5 GeV the 
size of A can exceed unity and the m h = 300 GeV curve approaches 3 at m D ~ 3 GeV. This 
may seem to signal the breakdown of perturbativity in the low-m D range, but an investigation 
into the perturbative unitarity of darkon-Higgs interactions at tree level JpJ has come up with 
the limit |A| < An ~ 12.6. Furthermore, it has been suggested in Ref. [7| that, due to a lack 
of clear division between the perturbative and nonperturbative regions of the parameter space, 
a reasonable requirement is |A| < 2a/7t (m^/100 GeV) 2 , which is roughly comparable to the 
preceding limit for the Higgs masses we have picked. Fourth, although the A values tend to 
become small as m D enters the region between 50 and 200 GeV or so, they get large again, 
approximately linearly with m D , as m D grows sufficiently large. This follows from the facts 
that (o" ann f rel ) is roughlv constant for the m D range of interest and that o" ann f rel ~ A 2 / {Anm 2 -^ 
for m D S> m-wzht' U* cL Il0 |- Lastly, the curves in the SM4+D appear quite similar to the 



corresponding ones in the SM3+D, but at most of the m D values considered we find the former 
to be lower than the latter. This decrease is mainly less than 20%, but it reaches nearly 25% 
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FIG. 1: Darkon-Higgs coupling A as a function of darkon mass m D for Higgs mass values m h = 
115, 200, 300 GeV in (a) SM3+D and (b) SM4+D with m t , = 500 GeV. The band widths in all figures 
correspond to the relic-density range which we have taken, 0.1065 < n D h 2 < 0.1181. 
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at m D ~ 60 GeV. The reason for the decrease is that the Higgs total width in the SM4 is, as 
mentioned earlier, enlarged relative to that in the SM3, which is also true for the total width 
of h in Eq. fll]). The enlargement ranges mostly from a few percent to ~40% and gets as high 
as ~70% at m h = 2m n ~ 120 GeV. 



III. CONSTRAINTS FROM DARK-MATTER DIRECT SEARCHES 



A number of underground experiments have been and are being performed to detect DM 
directly by looking for the recoil energy of nuclei caused by the elastic scattering of a WIMP 



off a nucleon 27H32|. Although indirect DM searches have recently turned up some intriguing 



results which may be interpreted as evidence for WIMPs 33j, it is very difficult to establish 
a firm connection to DM due to the indirect nature of the observed events. Therefore, direct 
detection is crucial to determine the properties of DM. 

In the SM4+D, the WIMP-nucleon interaction occurs via the exchange of a Higgs boson 
between the darkon and the nucleon N in the t-channel process DN — > DN. Thus to evaluate 
this elastic scattering requires knowing not only the darkon-Higgs coupling A, but also the 
Higgs-nucleon coupling g NNh , which parametrizes the Higgs-nucleon interaction described by 
£-NNh = ~9NNh-NN h. From this Lagrangian and C D in Eq. (j3J), one can derive for \t\ <C m\ 
the darkon-nucleon elastic cross-section [JO, d, @] 



* 9nnh v m N 



(5) 



7T 



(m D + m N ) m\ 

having used the approximation (p D + p N ) 2 — (m D + m N ) 2 . 

It remains to determine the value of g NNh , which is related to the underlying Higgs-quark 
interaction described by C qqh = —Yl q m q qq h/v, where in the SM4 the sum runs over the eight 
quark flavors, q = u,d, s,c,b,t,b',t' . Since the energy transferred in the darkon-nucleon scat- 
tering is very small, of order tens of keV, one can employ a chiral-Lagr ang ian approach to 
estimate g NNh - This has been done previously in the context of the SM3 [9|, [34j, [35| . In the SM4 
case, we modify the derivation described in Ref . j^] , incorporating the t' and b' contributions, to 
arrive at 



„SM4 
tJNNH 



III 



N 



17 

27 m B 



(6) 



where m N is the nucleon mass and m B denotes the lightest octet-baryon mass in the chiral limit, 
which can be related to the pion-nucleon sigma term, a nN , by m B ~ —13.39 cr^jy + 1.269 GeV 0. 
With a. 



nN 



45 MeV 36 , we obtain 



_SM4 
!JNNh 



2.11 x 10~ 3 , (7) 

to be compared with the SM3 value g^^ h = 1.71 x 10~ 3 0. We adopt these numbers in our 
numerical calculations below. We note, however, that a nN is not well determined, with values 
ranging roughly from 35 MeV to 80 MeV having been quoted in the literature 35M37||. which 
translate into 1.8 x 10~ 3 < gf^ A h < 3.3 x 10~ 3 and 1.3 x 10~ 3 < gfw h < 3.2 x 10" 3 . 
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With A and 

9nnh known, we can now predict the darkon-nucleon elastic cross-section cx el as 
a function of darkon mass once the Higgs mass is specified. We show our results for a el in Fig. [2], 
where the choices of darkon and Higgs masses are the same as those in Fig.[TJ For comparison, 
we display a el graphs in the SM3+D and in the SM4+D with m t , = 500 GeV. The SM4+D 
plots for m t , = 400 and 600 GeV are again similar to the one displayed. We find that at most 
of the m D values considered the a el curves in the SM4+D for a fixed m h are higher than the 
corresponding ones in the SM3+D. This difference ranges from a few to roughly 50 percent and 
results from the gf^ h enhancement relative to g^ff h overcoming the Asm4+d suppression relative 
to Asm3+d mentioned in the preceding section. 
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FIG. 2: Darkon-nucleon elastic cross-section a el as a function of darkon mass m D for Higgs mass values 
m h = 115, 200, 300 GeV in (a) SM3+D and (b) SM4+D with m t , = 500 GeV, compared to 90%-C.L. 
upper limits from XENON10 (black dashed-curve) , CDMSII (brown [gray] dashed-curve) , and CoGeNT 
(purple dashed-curve), as well as projected sensitivities of SuperCDMS at Soudan (green dot-dashed 
curve), SuperCDMS at Snolab (brown [gray] dotted curve), and XENON100 (black dotted curve). 
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In Fig.[2]we also plot the 90%-C.L. upper-limit curves for the WIMP-nucleon spin-independent 



elastic cross-section reported by the XENON10 28], CDMSII 30], and CoGeNT mjexper- 



Laep 

3L 



iments, along with the expected sensitivities of a number of future experiments [38| . For 



m D <10GeV, there are also limits from CRESST-I [27J and TEXONO [29j, but they are both 
above the predictions of the two models. 

Comparing the prediction curves of both models to the experimental upper-bounds in Fig. [2j 
one can see that some portions of the darkon mass regions considered are excluded, but the 
greater part of them are still viable. More precisely, for m h = 115, 200, and 300 GeV the 
XENON10 and CDMS II limits have ruled out darkon masses from ~9 GeV to between 70 and 
80 GeV, except for the 50 GeV < m D < 70 GeV range in the m h = 115 GeV case. Moreover, 
in the low-m D sections of the plots the exclusion limit from CoGeNT can be seen to rule out 
part of the 4GeV<m D <5GeV range. In contrast, darkon masses larger than 80 GeV or so 
are not yet probed by the current data from direct searches. As the projected sensitivities of 
future experiments in this figure suggest, SuperCDMS at Snolab and XENON100 may probe 
these models further to m D ~ 400 GeV, but SuperCDMS at Soudan may be unlikely to provide 
much stronger constraints on the models than the present bounds. 

It is interesting to point out that these two darkon models can accommodate the possibility 
that the excess events observed by CoGeNT originate from interactions with a relatively light 



WIMP of mass between 7 and 11 GeV [3jj, which is compatible with the two signal-like events 
detected by CDMS II [30j if they are also interpreted as evidence for WIMP interactions. For 
m D values within this range, the prediction curves in Fig. [2] each have some overlap with the 
possible signal region reported by CoGeNT [3l|. This is more so if we take into account the 
uncertainties in g NNh noted above, which could imply an increase in the predicted a cl by up to 
a factor of 3. 

Before moving on, it is worth remarking that, as Fig. [2] indicates, a el for a fixed m h approaches 
a constant value as m D becomes much greater than m wz h t ,. The reason is that in this large-m^ 
limit the ratio A 2 /m 2 D is, as pointed out in the preceding section, approximately constant and 
cr el in Eq. ([5]) is proportional to the same ratio, X 2 /m 2 D . Another observation from this figure is 
that the asymptotic value of <r el decreases as m h increases, which is in accord with Eq. ([5]). Hence 
direct DM searches in the future may lack the sensitivity to probe the larger darkon masses if 
the Higgs mass is also large. 



IV. SOME IMPLICATIONS FOR HIGGS SEARCHES AT COLLIDERS 

Since both the SM3+D and SM4+D have only a small number of free parameters, the relevant 
ones here being A, m D , and m h , it is possible to draw strong correlations among them [lol ]. This 
implies that these darkon models have a high degree of predictivity and that there is some 
simplification in testing them, without requiring many different observables. To illustrate this, 
we now discuss the Higgs decay into a pair of darkons, h — > DD, and some of its consequences 
for Higgs studies at colliders, in light of the bounds obtained above from comparing with DM 
direct-detection data. 
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Using the A values obtained in Sec.lIIIt we compute the rate and branching ratio of the 
invisible mode h — > DD. The results are depicted in Fig.[3l where the Higgs and darkon mass 
choices are the same as those in Fig.[TJ One can observe that the situations in the SM3+D and 
SM4+D are similar, namely that the sizable values of A in Fig. [1] translate into huge enhancement 
of the Higgs width via the additional process h — > DD and, consequently, an invisible branching 
ratio that is large. This is especially so if 2m D < m h < 2m w , in which case the Higgs partial 
width into standard particles is small. Although the 4th-generation quarks can cause the decay 
mode into a gluon pair, h — > gg, to dominate if m h < 140 GeV [15] , the inclusion of the darkon 
in the SM4+D leads to the dominance of h — > DD instead. 

The potential importance of the Higgs invisible decay mode in the darkon presence can 
alleviate some of the restrictions on Higgs masses in the SM4. For instance, at hadron colliders 
the important channel gg — >• h — )■ WW — > iviv is expected to be enhanced due to the new 
quarks in the SM4 by a factor of 9 for 100 GeV < m h < 200 GeV, the measurement of which would 



also provide indirect evidence for the new quarks 13J. Preliminary searches at the Tevatron for 



this channel have so far come back negative, thus excluding a large portion of this m h range 



in the SM4 13|, |39J. In the SM4+D, however, the possible dominance of the h — > DD mode 
implies that the enhancement of gg — >■ h — > WW — > iviv would likely be reduced or even 
negated completely, and so the Higgs-mass constraints could be weakened or evaded. 




50 

m D (GeV) 



50 
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FIG. 3: Partial width and branching ratio of invisible decay h —> DD as functions of darkon mass m D 
for Higgs mass values m h = 115, 200, 300 GeV in (a,c) SM3+D and (b,d) SM4+D with m t , = 500 GeV. 
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As discussed previously [10[ , collider measurements of the Higgs invisible decay with sufficient 
precision can lift some possible ambiguities in determining the darkon mass from direct DM 
searches. A substantial Higgs invisible decay can also be advantageous for testing the darkon 
models if the constraints from direct DM searches are combined with Higgs studies at colliders. 
As found above, the greater part of the darkon mass range from about 9 to 80 GeV in the SM3+D 
and SM4+D have been ruled out by direct-detection data if the Higgs mass m h = 115, 200, or 
300 GeV, the main exception being the neighborhood of m D ~ 57.5 GeV. Since Fig.[3] shows 
that a Higgs boson with one of these masses decays dominantly or significantly into a darkon pair 
if m D < 50 GeV, then the observation of such a Higgs boson with a sizable invisible branching- 
ratio might hint at inconsistencies of the models. All this illustrates that the interplay between 
direct DM searches and the study of the Higgs boson at colliders can yield crucial information 
about the darkon properties. 

An enhanced Bih — > DD) affects not only collider searches for the Higgs boson, but also 
decays which are mediated by or produce it. In the following two sections, we deal with such 
processes arising from the Higgs flavor-changing neutral couplings to quarks. 



V. CONSTRAINTS FROM B -> KDD DECAYS 



As seen in Sec. lIIIt direct DM searches with underground detectors currently being done or 
to be done in the near future are not expected to be sensitive to darkon masses below a few 
GeV. It turns out that such darkon masses can be probed using the decays of mesons containing 
the h quark. In this section we explore constraints available from the 5-meson decay B — > KDD, 
which contributes to the B decay into K plus missing energy, B — > K$. One could carry out 
a similar analysis using B — > K*DD, but we will not do so here. We will also briefly comment 
on the spin-one bottomonium decay T — > DD~f. 

Since the Higgs boson h is the only SM particle to which D couples, B — > KDD is induced 
by the flavor- changing 6-quark decay b — > sh* — > sDD, the effective bsh coupling being loop- 
generated with up-type quarks and the W boson in the loops. These transitions have been 
studied previously in the context of the SM3+D in Refs. 0, Generalizing their results to 
the SM4+D, we can express the effective Hamiltonian for b — > sh* — > sDD as 

= ^f^l^)^ 2 , (8) 

where 

a* = £W S ^ + A ^') , x q = £| , x b q s = v; s v qb , (9) 

with Vki being the elements of the 4x4 Cabibbo-Kobayashi-Maskawa (CKM4) matrix, and con- 
tributions from u and c quarks have been neglected. Hence the corresponding expression for g bs 
in the SM3+D does not contain the \ b t fx t , term. 
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The amplitude for B — > K DD is then 



M(D-^h-DD) = *^»/ (S), (10) 



where fn(s) = 0.3 exp(0.63 s/m 2 B — 0.095 s 2 /m B + 0.591 s 3 /m^) is the relevant B — > K form- 
factor [2J, with s = (p B — Vk) 2 i an d the approximation (m B — m^A / (m b — mj ~ m B /m b has 
been made. It follows that 

r (B +^+ DD ) = -|^/K), (ii) 

where a factor of 1/2 has been included to account for the identical .D's in the final state and 



J(m D ) = / dS(/ (S)) V 1 -— -VK"^" j ) -4m^5. (12) 

J4m 2 D ' s 

We note that our formula for the I? — )■ KDD rate agrees with the corresponding one obtained 
in Ref. [12}, but is 4 times smaller than that given in Ref. Q. 2 Since for m D <C m ft we can 
simplify Eq. @ to 

°ann u rel — 4 ™ ' V ±0 / 

ml mo 



we can rewrite Eq. ( fill as 

1 ; 2048vr 3 t; 2 ^(h^Xi) ' 1 ' 

where both cr ann t> rel and T,jT(h — » Xj) have m D dependence. We can get constraints on the 
parameter space m D < (m B — m K )/2 by comparing this prediction with the experimental 
information on the B decay into a kaon plus missing energy, which receives contributions from 
B — > KDD and B — > Kvv. We first update the constraints in the SM3+D and then discuss 
the SM4+D case. 

The prediction of the branching ratio B(B — > KDD) in the SM3+D using Eq. (Till) , with the 
\\fx t , term in g bs dropped, involves large uncertainties which come mainly from the calculation 
of the total width Y,iT(h — > X^j for m~ h under a few GeV. In the case of the physical h, for 
m h < 2 GeV the predicted rate of the important channel h — > hadrons is well known to contain 



significant uncertainties [40j, |4l|. In estimating the h total width, for 2m n < m~ h < 1.4 GeV 



we adopt the T(h — > hadrons) results from Ref. 4JJ, whereas for smaller and larger values 



of m~ h we simply use the perturbative formulas for Higgs decays 42|. We graph the resulting 



B(B — > KDD) as a function of m D in Fig. HI which is to be compared with experimental data. 



2 This can be traced to a factor of 1/2 apparently missing in the expression for the B — > KDD amplitude in 
the Eq. (6) of the first paper in Ref. 
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FIG. 4: Branching ratio of B + — > K+DD as a function of darkon mass m D in SM3+D 



The latest experimental search for the decay B+ — > K+vv has produced the branching-ratio 
limit B cxp (B+ — > K+vv) < 14 x 10~ 6 [43]. On the theoretical side, the most recent calculations 
in the SM3 predict B(B+ — > K+vv) to be between 3.6 x 10 -6 and 5.1 x 10~ 6 , with errors of 
order 15% [441 ] . Accordingly, it is reasonable to take B exp (B+ — > K+~$) ~ B exp (B+ — > K+vv), 
from which we can subtract the SM3 prediction for B(B+ — > K+vv) in order to require that 
B[B+ K+DD) < 1 x 10~ 5 in the SM3+D. 

We can see from Fig.Hthat the predicted B(B+ — > K+DD) > 1 x 10~ 5 for m D < 1.5 GeV. 
Recalling the hadronic uncertainties mentioned above, we can then conclude that in the SM3+D 
much of this range of m D values, especially m D < 0.4 GeV, is excluded by the data. We would 
need improved data from future measurements of B — > Klfi before we could disallow more 
darkon masses within the m D < 2.4 GeV region. These conclusions are similar to those made in 
Ref. [z| due partly to the stronger experimental limit at present and partly to the overestimate 
of their B(B ->• KDD). 

In the SM4+D, the prediction for T(B — > KDD) is modified due to the presence of the 
new quarks, t' and b'. The loop- induced effective coupling g bs in Eq. ffT4"|) receives a t'-quark 
contribution as given in Eq. ([9]). To examine its effect on T(B — > KDD), we need to compare 



„SM4 f n ^SMS 
ybs LO iJbs 



191 from 



sM4/ fl SM3|2 for m</ = 4 oo a nd 500 GeV 
600 GeV. In addition, the presence 



For concreteness, we take the relevant CKM4 elements extracted in Ref. 
a global fit for the SM4. Accordingly, we can expect that the numbers we use are typical values 
for the model. Thus, with \\ s = 0.04 in gf s M3 , we find \gf 
to be similar in value, ~1.2, but it goes up to 1.6 for m t , 
of t' and b' affects the total width SjF(/i — > X^j in Eq. ( TT4l) mainly via their contributions to 
h — > hadrons due to the quark-loop induced h — > gg, as already mentioned earlier. Despite 
the hadronic uncertainties, this implies that the enhancement of the h total width in the SM4 
compared to the SM3 can be expected to be less than (5/3) 2 ~ 2.8 if 2m n < m~ h < 2m c , where 
5 and 3 are the numbers of heavy quarks in the two models, respectively, for this m~ h range. This 
enhancement decreases to no more than 25% after the h — )■ cc channel is open. We can then 
conclude that the effects of these two factors on T(B — > KDD) in the SM4+D amount to changes 
to the rate in the SM3+D by no more than a factor of 2 in either direction, implying that the curve 
in Fig. H] would not be very different in the SM4+D. Since the prediction for B(B + — > K+vv) is 
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raised by at most 20% in the SM4 [lSfl , the empirical bound on B(B -> KDD) in the SM4+D 
can be taken to be unchanged compared to that in the SM3+D. It follows that the constraints 
on the darkon masses within the m D < 2.4 GeV range in the SM4+D are similar to those in 
the SM3+D. 

We should also mention that for m D <170MeV an additional restriction is provided by 
the kaon decay K + — > ti + $, which receives a contribution from K + — > tt + DD. The agree- 
ment between the SM3 expectation and experimental data on K + — > implies that for 
m^, < 170 MeV the SM3+D is already ruled out pj, as is the SM4+D. This is consistent with 
what can be inferred from the low-m D end of Fig. HI 

For the larger range m D < 5 GeV, there may also be constraints available from future measure- 
ments of the decays T — > DD^y. Presently, the existing experimental limits on T — > 7 + invisible 



are not yet strong enough to probe these darkon models 11 



VI. FCNC DECAYS Q qDD 

The presence of the new quarks in the SM4+D can have important implications for probing the 
darkon sector that are lacking or absent in the SM3+D. In the SM3 the flavor-changing neutral 



45, 


46|, 


19, 


47 



We expect that in the SM4+D the 
related decay t — > ch* — > cDD, if kinematically allowed, can be similarly enhanced. These 
processes may be detectable at the LHC after its operation reaches full capacity in the near 
future. The Tevatron and the LHC can also produce the new quarks, t' and b', if they exist, in 
a similar way as they can produce the t quark, albeit fewer of them due to their bigger masses. 
It is therefore of interest as well to explore their FCNC decays, t' — > (c,t)h* — > (c,t)DD and 
b' — > (s,b)h* — > (s,b)DD, which may have observable rates. These decays could, in principle, 
probe darkon masses from zero all the way up to (itlq — m q )/2, hence covering potentially wider 
m D ranges than those covered by some of the DM direct searches in the future. Here we estimate 
the branching ratios of these FCNC decays involving the darkon. The corresponding decays with 
the u and d quarks, — > uDD and b' — > dDD, are comparatively suppressed due to the less 
favorable CKM4 factors. 

The Lagrangian describing the FCNC transition Q — )■ qh involving a heavy quark Q and 
a lighter quark q can be written as 

C Qqh = q(g? q P L + gl q Pn)Qh, (15) 

where P LR = |(1 =F 7 5 ) and the loop-induced couplings g® q R generally depend not only on the 
internal quark (and W) masses and the CKM matrix elements, but also on the masses and 
momenta of the external particles. The amplitude for Q — > qh* — > qDD is then 

M(Q -> qDD) = n , L , 16 

m z h -s- iY h m h 
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where s = (p Q - p q 
r(Q ->• gDD) 



This yields the decay rate 



2567T 3 77lg 



(m -m„) 2 



1 - 



4m 



m? — s)' 



x 



(m 2 Q + m\ 



4m 2 , s 



s)+4Re(g^g^)m Q m q 



(s-m 2 h ) +T 2 h 



(17) 



m 



For 2m D < m h < rriq — m , the Higgs-pole contribution dominates this integral, and so one has 
r(Q -»■ gDD) ~ r(Q -»■ g/i) -D-D). 

The effective couplings g® ^ in the SM have been evaluated previously for arbitrary values of 

We make use of the formulas provided in Ref. 49 



the external and internal masses 47-49 



In our numerical illustration below, we take m h = 115 GeV and m# = 500 GeV, as well as 
the corresponding elements of the CKM4 matrix extracted from a global fit in Ref. 19] . To 
determine the branching ratios, we normalize the decay rates according to 



B(t cDD) 
B{t' -> qDD) 
B(b' -> gDD) 



T(t -> cD£>) 

Tit^bwY ' 

r(f -»■ gDD) 
T(t' -+ bW) + r(t' sW) 
r(6' ->■ gDD) 



(19) 
(20) 



r(6' -> + T{V -> cW) ' 

following Ref. 19[ in the Q qh cases. We display the results in Fig. |5j 

Estimates suggest that t — > ch can be detected at the LHC if its branching ratio is several 
times 10~ 5 or higher [46}. In the presence of the darkon, if h — > DD is leading, t — > ch — > cDD 
is more likely to occur than other t — > ch — > cX modes. In that case, however, as Fig. [5] 
indicates, B(t — > cDD) < 1.0 x 10~ 8 and so t — > cDD probably will not be observable in 
the near future. The branching ratio could be several times higher if m t > ~ 700 GeV, but this 



already exceeds the perturbative unitarity upper-bound m t i ~ 550 GeV [13 



0.1 
0.001 
f 10" 5 
lO" 7 
10"" 



1 



t'^tDD 



t'^cDD 



t^cDD 



(a) 
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FIG. 5: Branching ratios of (a) t -> cDD, t' -> cDD, and t' -> tDD and (b) b' -> sDD and b' -> bDD 
as functions of darkon mass m D for Higgs mass m h = 115 GeV in SM4+D with m t , = 500 GeV. 
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In contrast, the t' and b' numbers in Fig. [5] are much greater: B(t' — > cDD) < 8.2 x 10 5 , 
B(f -»■ WD) < 1.4 x 10~ 3 , £(&' -»■ sDf ) < 3.1 x 1(T 4 , and B(b' bDD) < 3.3 x 10~ 3 . Since 
t' — > qh and b' — > qh decays with branching ratios between 10~ 4 and 10~ 2 are expected to be 



within the reach of the LHC |19|, |47j, we may expect that these if — > qDD and b' — > qDD 
decays would also be detectable at the LHC despite their final states involving missing energy. 
Once they are measured, comparing the results with those from DM direct searches could provide 
additional consistency tests for the darkon models. 



VII. CONCLUSIONS 



We have explored one of the simplest dark-matter models, the SM4+D, consisting of the 
standard model with four generations and a real gauge-singlet scalar, the darkon, to play the 
role of WIMP dark matter. This model possesses not only the phenomenologically interesting 
features of the SM4, but also a high degree of predictivity in its DM sector. We have investigated 
constraints on the SM4+D from DM direct-search experiments and from S-meson decay into 
a kaon plus missing energy. Compared to the SM3+D case, the resulting bounds in the SM4+D 
are similar, namely that for the representative Higgs masses chosen most of the darkon masses 
between roughly 4 to 80 GeV are excluded by the direct searches and that much of the mass 
region below 1.5 GeV is also excluded by the B decay data. Interestingly, the SM4+D as well 
as the SM3+D can also accommodate the possible interpretation that the excess events recently 
measured by the CDMS II and CoGeNT experiments were due to interactions with a light WIMP 
of mass around 9 GeV. Darkon masses greater than 80 GeV in the two models are still viable 
and can be probed by future direct searches. 

We have discussed the complementarity of DM direct searches and Higgs studies at colliders 
in testing the darkon sector of the SM4+D. This can be crucial for a relatively light Higgs boson, 
which may decay substantially into the invisible darkons. Accordingly, we have pointed out that 
existence of the darkon could lead to the weakening or evasion of some of the restrictions on the 
Higgs mass in the presence of fourth-generation fermions. 

We have considered some implications of the SM4+D that are lacking or absent in the SM3+D 
as far as probing the darkon properties is concerned. In particular, we have examined the 
Higgs-mediated FCNC decays t — > cDD, t' — > (c,t)DD, and b' — > (s,b)DD, which may 
have observable rates at current or future colliders. These processes promptly proceed from 
the Q —7- qh transitions if the decay mode h — > DD is dominant. Although the t — > cDD 
branching-ratio is enhanced by several orders of magnitude compared to that in the SM3+D, 
reaching the 10 -8 level, this decay is still unlikely to be measurable in the near future. In contrast, 
the branching ratios of t' — > qDD and b' — > qDD can be as large as a few times 10~ 3 , which 
may be detectable at the LHC. If observed, they would offer extra means to test the models, 
covering darkon masses from zero up to hundreds of GeV. 
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